(N 
O 



Selection rules in the (3(3 decay of deformed 

nuclei 



S-H 
OH 

^ ■ Jorge G. Hirsch * Octavio Castanos } Peter O. Hess * 

Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
A. P. 70-543 Mexico 04510 D.F. 

Osvaldo Civitarese § 

Departamento de Fisica, Universidad Nacional de La Plata, 



o 
o 
in 
o 



' c.c.67; 1900, La Plata, Argentina 



O 
=5 



February 8, 2008 



Abstract 



The 2u (3(3 decay half-lives of six nuclei, whose decays were previ- 
■ ously reported as theoretically forbidden, are calculated by including 

the pairing interaction, which mixes different occupations and opens 
up the possibility of the decay. All allowed channels for the 0v (3(3 de- 
cay are also computed. The estimated 2u (3(3 half-lives suggest that 
measurements in 244 Pu may find positive signals, and that planned ex- 
periments would succeed in detecting the Pfou decay in 160 Gd. Limits 
for the zero neutrino mode, in the analyzed deformed emitters, are 
predicted. 

Pacs numbers: 21.60.Pw, 23.40.Hc, 27.70. +q, 27.90. +b 

Neutrinoless double beta decay (/3/3o„), if detected, would offer definitive 
evidence that the neutrino is a Majorana particle, i.e. that it is its own 
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antiparticle |], |2|. It would also provide the information needed to deter- 
mine neutrino masses, complementary to the one obtained from solar and 
atmospheric neutrino experiments ||. £| . 

Theoretical nuclear matrix elements are needed to convert experimental 
half-life limits, which are available for many /5/5-unstable isotopes 0, [J, into 
constrains for the effective Majorana mass of the neutrino and the contribu- 
tion of right-handed currents to the weak interactions. Thus, these matrix 
elements are essential to understand the underlying physics [||, [7], The 
two neutrino mode of the double beta decay {fifiiv) is allowed as a second 
order process in the standard model. It has been detected in ten nuclei || 
and it has served to test a variety of nuclear models || . 

The pseudo SU(3) approach has been used to describe many low-lying 
rotational bands, as well as B(E2) and B(M1) intensities, in rare earth and 
actinide nuclei, both with even-even and odd-mass numbers. The theoretical 
results show, in general, a very good agreement with the data 0. The (3(3 
half-lives of some heavy deformed nuclei, which may decay to the ground and 
excited states of the daughters, were evaluated for the two and zero neutrino 



emitting modes |I| [TTJ, [L2|, [13], |14| using the pseudo SU(3) scheme. The 



predictions were in good agreement with the available experimental data for 
150 Nd and 238 U. The double electron capture decay channel was studied for 
the decay of other three nuclei |L5| . 

The simplest pseudo SU(3) model predicts the complete suppression of 
the (3(3 2v decay for the following five nuclei: 154 Sm, 160 Gd, 176 Yb, 232 Th and 
244 Pu |Ll], [14]]. Recently, it was argued that the cancellation of the (3(3 2v decay 
in 160 Gd would suppress the background for the detection of the Ou mode 



In the present contribution we extend the previous research [1C, 11, 12 



g |TJ, and evaluate the (3(3 half-lives of 154 Sm, 160 Gd, 170 Er, 176 Yb, 232 Th 
and 244 Pu using the pseudo SU(3) model. In these nuclei the 2v(3(3 mode is 
forbidden when the most probable occupations are considered. To be able to 
evaluate finite half-lives, we were forced to include in the calculations states 
with different occupation numbers which can be mixed through the pairing 
interaction. The amount of this mixing is evaluated, and the possibility of 
observing the (3(3 decay is discussed for both the 2v and Ov modes. 

Pseudo-spin symmetry |T8|] describes the quasi- degeneracy of the single- 



particle orbitals with j = 1 — 1/2 and j — (I — 2) + 1/2 in the rj shell. It 
allows these orbitals to be classified as pseudo-spin partners with quantum 
numbers j=j,fj = i] — 1 and 1 = 1 — 1. The first step in the pseudo SU(3) 
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Table 1: Deformations, proton and neutron occupation numbers, pairing mix- 
ing hp a i r and excitation AE energies (both in MeV) and mixing coefficients 

description of any nuclei is to find the occupation numbers for protons (p) 
and neutrons (n) in the normal and abnormal parity states rip,n^,n A ,n A . 
These numbers are determined by filling the Nilsson levels from below, as 
discussed in [[HJ. The deformations (19| and occupancies for the 12 isotopes 



studied in the present work are shown in Table 1. 

In the f3(32 V decay each Gamow- Teller operator annihilates a proton and 
creates a neutron in the same oscillator shell and with the same orbital 
angular momentum. As a consequence, the fifav decay is allowed only if the 
occupation numbers fulfill the following relationships 

n p,f = n p,i + ^ J n n,f = n n,i ' 

n p,f = n p,i J n nJ = n n,i — 2 . (1) 



This selection rule [10] forbids the fifiiv decay between the nuclei marked 
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with (1) or without comments in Table 1. However, the pairing interaction 
allows the mixing between states in the same nuclei with pairs of nucleons 
transferred between different configurations. These excited configurations are 
indicated by (2) in Table 1. An energy difference AE is required to promote 
a pair of nucleons from the last occupied normal parity orbital to the next 
intruder orbital (or viceversa), in the deformed single particle Nilsson scheme. 
It is listed in the seventh column of Table 1. 

The "leading SU(3) irreps", those which are the most bounded under 
the quadrupole - quadrupole interaction, are the dominant component of the 
ground state wave function in these heavy deformed nuclei, representing in 
most cases more than 60% of the total wave function ||. For the eight rare 
earth isotopes listed in Table 1, this dominant wave function component can 
be written as 



The proton, neutron, and total SU(3) irreps associated to each set of 
occupation numbers are listed in Table 2. 

As a first approximation, we will describe the ground state of each nucleus 
CIS db linear combination of these two states: 



with x\ + x\ = 1. 

Many multipole-multipole pairing type interactions can remove a pair of 
nucleons from an unique parity orbital and create another pair in a normal 
parity one. In the present approach we are restricting the pairs of nucleons in 
intruder orbits to be coupled to J=0, i.e. to have seniority zero. Under this 
approximation the only term in the Hamiltonian which can connect states 
with different occupation numbers in the normal and unique parity sectors 
is pairing. In the present case, the Hamiltonian matrix has the simple form 



with h pair = (Nucleus(2), + |f/p a j r |Nucleus(l), + ), whose explicit expression 
is given in |17|]. The values of h pa i r ,xi and x<i are presented in the last two 




Nucleus, + ) = xi |Nucleus(l), + ) + x 2 |Nucleus(2), + ) 




(3) 
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Nucleus 


(Ap,//p) (A„,//„) (\,n) 


154 Sm (1) 
154 Sm (2) 
154 Gd 


(10, 4) (18, 0) (28, 4) 
(10, 4) (18, 4) (28, 8) 
(10, 4) (18, 0) (28, 4) 


ibo Gd 

160 Dy (1) 
160 Dy (2) 


(10, 4) (18, 4) (28, 8) 
(10, 4) (18, 0) (28, 4) 
(10, 4) (18, 0) (28, 4) 


i7u Er 
170 Er (2) 

170 Yb 


(10,4) (24, 0) (34, 4) 
(4,10) (24, 0) (28, 10) 
(4,10) (20, 4) (24, 14) 


176 Yb 

176 Hf (1) 
176 Hf (2) 


(4, 10) (18, 8) (22, 18) 
(0, 12) (20, 6) (20, 18) 
(4, 10) (20, 6) (24, 16) 


232 Th (1) 
232 Th (2) 

232 jj 


(12, 2) (30, 4) (42, 6) 
(18, 0) (30, 4) (48, 4) 
(18, 0) (26, 4) (44, 4) 


244 Pu (1) 
244 Pu (2) 
244 Cm 


(18, 0) (34, 8) (52, 8) 
(18, 4) (34, 8) (52, 12) 
(18, 4) (34, 6) (52, 10) 



Table 2: The proton, neutron, and total irreps assigned to each nucleus 
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columns of Table 1. In the case of 244 Pu we are using a small deformation 



[ 19| . for which the two configurations listed in Table 1 are nearly degenerate, 
and have maximal mixing. It has important consequences upon its predicted 
(3(3 half-life. 

The inverse half- life of the two neutrino mode of the /3/3-decay, /3/3 2v , can 
be evaluated as [EO 



4 /2 (o- 



H 



G 2v I M, 



2v 



(4) 



where G 2 „ is a kinematical factor which depends on Qpp, the total kinetic 
energy available in the decay. 
The nuclear matrix element is 



M- 



2v 
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GT 
2v 



(o+||r||i+) (i&||r||o+) 



N 



E f + E. 



N 



Eh 



(5) 



being T the Gamow- Teller operator. The energy denominator contains the 
intermediate E N , initial Ei and final Ef energies. The kets denote 
intermediate states. 

The mathematical expressions needed to evaluate the nuclear matrix el- 
ements of the allowed g.s. — > g.s. (3(3 decay in the pseudo SU(3) model were 
developed in flO|. Using the summation method described in JTO 



exploiting 

the fact that the two body terms of the SU(3) Hamiltonian commutes with 
the Gamow- Teller operator [11 1, resuming the infinite series and recoupling 



the Gamow- Teller operators, the following expression was found: 



Mi 



GT 
2u 



°jp, n) 
£ 



[al ® an] 1 



[al®a n ] 1 



j=o 



(6) 



where a(p, n) are the 1-body Gamow- Teller matrix elements and the energy 
denominator £ is determined by demanding that the excitation energy of the 
Isobaric Analog State equals the difference in Coulomb energies [ 10|, [17| . The 



index p, n refer to the orbitals %/ 2 ,«i3/2 f° r Sm, Gd, Er, and Yb, 
and to the orbitals jf 3//2 , j™ 5 / 2 for 232 Th and 244 Pu. 

As it was discussed in fl(| Eq. (|]) has no free parameters, being the 
denominator £ a well defined quantity. The reduction to only one term 
comes as a consequence of the restricted proton and neutron spaces of the 
model. The initial and final ground states are strongly correlated with a very 
rich structure in terms of their shell model components. The evaluation of 
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the matrix elements in the normal space of Eq. (HD is performed by using 
SU(3) Racah calculus to decouple the proton and neutron normal irreps, and 
expanding the annihilation operators in their SU(3) tensorial components. 

For the six potential (3(3 emitters listed in Tables I and II, the flfov decay 
can only proceed through the second component of the ground state wave 
function, and for this reason it is proportional to the amplitude x 2 . Its 
explicit expression is given in [|17[ . 



For massive Majorana neutrinos one can perform the integration over the 
four-momentum of the exchanged particle and obtain a neutrino potential, 
which for a light neutrino (m u < 10 MeV) has the form 

TTI — . 2R f°° , sin(qr) 

H(r, E) = — dq—^- , 7 
pr Jo q + E 

where E is the average excitation energy of the intermediate odd-odd nucleus 
and the nuclear radius R has been added to make the neutrino potential 
dimensionless. In the zero neutrino case this closure approximation is well 
justified pi|| . The final formula, restricted to the term proportional to the 



neutrino mass, is |E| ^ 



(to?)- 1 = j GoX , (8) 

where G$ v is the phase space integral associated with the emission of the two 
electrons. The nuclear matrix elements M 0u are [2C] 



M , = \Mg T - 4i<J, with Mq v = (0+||Ol0+>, (9) 
9a 

where the kets |0^~) and |0j ) denote the corresponding initial and final nuclear 
states, the quantities gy and qa are the dimensionless coupling constants of 
the vector and axial vector nuclear currents, and 



gt = Y j B mt - m .a n t~H(\r m -r n \,E) , (10) 

m,n 

F = Y J t- m tnH{\r m -r n \ ) E) , 

m,n 

being a the Pauli matrices related with the spin operator and t~ the isospin 
lowering operator, which satisfies t~\n) = \p). The superscript GT denotes 
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154g m 


1.251 


4.872 10 


-21 


0.0445 


1.04 10 23 


160Q d 


1.730 


8.028 10 


-20 


0.0455 


6.02 10 21 


170 Er 


0.654 


6.496 10 


-23 


0.0374 


1.10 10 25 


176 Yb 


1.086 


3.866 10 


-21 


0.0306 


2.77 10 23 


232 Th 


0.858 


7.410 10 


-21 


0.0504 


5.30 10 22 


244 Pu 


1.352 


4.081 10 


-19 


0.0617 


6.43 10 20 



Table 3: The Q- values, phase-space integrals, matrix elements and predicted 
half-lives for the j3p 2 u beta decay 



the Gamow- Teller spin-isospin transfer channel, while F indicates the Fermi 
isospin one. In the present work we use the effective value {^) 2 = 1-0 f|. 

Transforming this operator to the pseudo SU(3) space, we arrive to the 
expression 

O a = 0% pNn + o% pAn + o ApNn + ApAn , (11) 

where the subscript index NN, NA, . . . are indicating the normal or abnor- 
mal spaces of the fermion creation and annihilation operators, respectively. 



In previous works |T(J |TT| we restricted our analysis to six potential dou- 
ble beta emitters which, within the approximations of the simplest pseudo 
SU(3) scheme, were also decaying via the 2v mode. They included the ob- 
served 150 Nd — > 150 Sm and 238 U — > 238 Pu decays. In these cases two neutrons 
belonging to a normal parity orbital decay in two protons belonging to an 
abnormal parity one. The transition is mediated by the operator ApNn . 

In Table 3 the (3 fa* decays of the six nuclei, previously reported as for- 
bidden, are presented. Those with the larger Qpp values have the larger 
phase-space integrals G 2v . The /^/^-decay matrix elements M^J are sup- 
pressed by a factor x 2 ~ 1/3 compared with the "allowed" ones (see [PUD, 
which reflects in ftfav half-lives a factor 10 larger than in other nuclei with 
similar Q- values. The exception is 244 Pu, which for the deformation used has 
a large mixing, and the shorter fifav half-life, which is not far from the limits 
reported in the Livermore experiment |p3 . The decay of 160 Gd is suppressed 



but it is still not far from the present limits [24], and large enough to be seen 
in the proposed experiments [16 |. 

Those configurations in which the /3/?2v transitions are forbidden can still 
be connected through the zero neutrino mode, due to presence of the neutrino 



8 



Nucleus 


G 0u [yr" 




M 0u 


TqL k m v) 2 [yr eV 2 ] 


154g m 


4.898 10" 


-15 


2.384 


9.38 10 24 


160 Gd 


1.480 10" 


-14 


0.919 


2.09 10 25 


170 Er 


1.673 10" 


-15 


0.731 


2.92 10 26 


176 Yb 


6.817 10- 


-15 


0.772 


6.43 10 25 


232 Th 


3.160 10- 


-14 


1.232 


5.44 10 24 


244 Pu 


1.463 10- 


-13 


1.171 


1.30 10 24 



Table 4: The phase-space integrals, matrix elements and predicted half-lives 
for the Oz/ double beta decay 



potential. In this way there are two terms in the PPo v decay: one connecting 
to the basis state which has allowed Pfov decay, and one to the state with 
forbidden ft feu decay. The equations needed in the first case are the same 



employed in the study of allowed decays [[10], [TTJ] . The second one involves 
the annihilation of two neutrons in normal parity orbitals, and the creation 
of two protons in normal parity orbitals (intruder- intruder in 154 Sm). The 
transition is mediated by the operator 0% Nn (0 ApAn ). A detailed description 
of the calculations involved is presented in fTTf . (3(3q u phase-space integrals, 
nuclear matrix elements and half-lives are shown in Table 4. 

As a consequence of the explicit inclusion of deformation in the present 
model, the f3(3o u half-lives are larger than those reported in [25]. In 160 Gd the 
/3/3o^-decay half-life is at least three orders of magnitude larger than the 
/?/3 2y -decay half-life. It implies that the background suppression due to a 
large /3/3 2 ^ half-life would be effective, although not as noticeably as was 



optimistically envisioned in ||16fl . In any case, the results presented strongly 
suggest that the planned experiments using GSO crystals [T(| would be able 
to detect the Pfov decay of 160 Gd, and to establish competitive limits to the 
(3(3q u decay. 

The present results consider only the dominant pseudo SU(3) irrep for 
each configuration. We have learned from realistic calculations, where the 
single particle term and pairing interactions induce the mixing of different 
irreps, that the leading irreps represent in most even-even heavy deformed 
nuclei at least 60 % of the total wave function ||. The inclusion of spin de- 
pendent terms in the Hamiltonian, relevant to the description of the Gamow- 
Teller resonance, is not expected to strongly modify the ground state wave 
function of the even-even initial and final nuclei. This dominance lead us 
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to expect that future calculations, which will take into account contribu- 
tions from various irreps, would slightly affect the present predictions. Given 
the leading role play by the quadrupole-quadrupole interaction in heavy de- 
formed nucleus, we are confident that the order of magnitude of the predicted 
(3/3 half-lives, when various irreps are included in the calculations, will remain 
unchanged, as compared to the results reported above. 

Work supported in part by CONACyT, and by a CONACyT-CONICET 
agreement under the project Algebraic methods in nuclear and subnuclear 
physics. O.Civitarese is a fellow of the CONICET, Argentina. 
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